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Abstract 

Green’s functions axe derived for the plane edaatoetatics problem of a dislocation 
in a bimaterial atrip. Using these fundamental solutions as kernels, various problems 
involving cracks in a bimaterial strip are analyzed using singular integral equations. For 
each problem considered, stress intensity factors are calculated fox several combinations 
of the parameters which describe loading, geometry and material mismatch. 

1 Introduction 

The mechanical behavior of bimaterial interfaces in composite materials is currently a topic 
of considerable interest to the applied mechanics community. Many analytical, numerical 
and experimental investigations have been conducted recently to gain a better understand- 
ing of bow these interfeces affect bulk composite properties such as strength, stiffness and 
toughness. The mechanisms of cracking or debonding along bimaterial interfaces are of 
particular interest in brittle composites since it has been established experimentally and 
through micromechanical models that certain desired properties can only be achieved in 
such materials if the cracks which initiate in the matrix are deflected by the fibers along the 
fiber-matrix interface [1]. It is clear that micromechanical analyses will continue to play an 
important role in analyzing and designing brittle composites to ensure this desired failure 


1 


sequence. Moreover, experimental programs are necessary to measure the fracture tough- 
ness of the fiber-matrix interface because the conditions required for the desired behavior 
involve the relative toughness between the fiber and the interface. 

Micromechanical models which involve relatively simple geometries such as cracks in 
infinite and semi-infinite plane bodies have been handled using distributed dislocations and 
singular integral equations [2]. For complicated finite geometries, on the other hand, the 
finite element method has gained popularity [3,4]. The singular integral equation method 
has two advantages. First, it leads to accurate results for stress intensity factors. Second, 
once the Green’s functions are derived and the equations are set up, parameter studies can 
be performed by simply varying the dimensionless parameters which describe the loading, 
crack length and geometry. Unfortunately, if the geometry of the problem is complicated, 
the method is not feasible because it is very difficult to derive the kernels. The real advantage 
of the finite element method is its ability to model complicated geometries. If proper care 
is taken, the method also produces accurate results. Parameter studies, on the other hand, 
are time consuming and relatively cumbersome. 

This paper addresses a class of problems which involve cracks in bonded strips. The 
motivation of this work came from a desire to develop a computer program which could 
be used to calculate stress intensity factoro and energy release rates for bimaterial fracture 
specimens being developed at C.W.&.U., and model the geometry of the so-called Santa 
Barbara mixed-mode specimen shown in Figure la. The latter “T-crack” configuration was 
developed by Charalambides et al. [3,4] at the University of California at Santa Barbara to 
measure the fracture resistance of bimaterial interfaces. To interpret experimental results 
and to guide the design of the specimen, they developed a finite element approach to charac- 
terize trends in stress intensity foctors, energy release rates and center paint displacements 
with specimen dimensions, elastic properties and crack length. As discussed in [3,4], the 
stress intensity factors and energy release rate for the crack growing along the interface 
exhibit steady state behavior as a result of the constant moment within the inner load- 
ing points. The shortest crack for which results were calculated in [3,4] is approximately 
a/l — 0.0938 (or a/k 2 = 0.3127). The results showed that this crack length is already 
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at steady state. One of the questions left unanswered is: how long does the crack have 
to grow along the interface before it reaches steady state? To answer this question the 
problem is modeled in this paper using the singular integral equation technique. Results 
will be presented for relatively short crack lengths as well as for the case a = 0, c < 
(Fig.lc, Fig.ld, Fig.le). It should be noted that at the end of this work the authors learned 
that Charalambides has recently obtained results for both the transient region and for the 
three-point loading configuration. 

In the next section Green’s functions are derived for edge dislocations in a bimaterial 
strip. These fundamental solutions can be used in turn to set up the integral equations 
for all the configurations shown in Figure 1. It should be noted that this class of problems 
has been analysed by Lu [5] and Lu and Erdogaa [6] using Fourier transforms. In [5,6] 
the problems were reduced to singular integral equations without the use of dislocations as 
fundamental solutions. Dislocation densities, however, were defined eventually to set up the 
integral equations. In this paper the dislocations are introduced as fundamental solutions. 
This will enable us, in the future, to solve problems which involve non-symmetric loading 
and/or inclined cracks. In the third section the dislocation solutions are used to set up 
the integral equations for several configurations and loadings. The last section presents 
numerical results and comparisons with existing solutions. 

2 Fundamental Solutions and Loading Conditions 

Assume that two dissimilar elastic strips are bonded along the x axis. The upper layer 
(y > 0) is labeled by “1” and the lower layer (y < 0) is labeled by “2”. The Green’s functions 
for an edge dislocation in the bimaterial strip are derived by superposing the solutions for 
(1) a dislocation near the interface of two bonded half-planes; and (2) a bimaterial strip 
loaded with boundary tractions which cancel out those induced by problem (1). 

The solution to problem (1) is well known. The stresses and displacements, which will 
be denoted by superscript “(l)”, can be expressed conveniently in terms of Muskhelishvili’s 
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complex potentials as 


(»w - i«W) I = »i + *,(2) + 1*5(2) + #i(2), 


(»a ) +^ ) ).=2[*iW+*iW]. 


Jw( iir +i isr )i= *♦<« - (•<« + *iw + *1 wi. 


(i) 


in which the subscript i ( i = 1,2) denotes "in region $1 and ¥1 correspond to the 
potentials for the npper half plane and $2 and $2 correspond to the potentials for the lower 
half plane. Moreover, z is the complex variable x + iy, the prime denotes differentiation 
with respect to z, an overbar denotes conjugation, /* is the shear modulus, and « is defined 
in terms of Poisson’s ratio v as * = 3 — \v for plane strain, and it = (3 - u)/( 1 + v) for 
plane stress. The complex potentials for a dislocation located at z& = a?o + »yo are given by 

(71 


$, = *? + #f , = *? + «? 


( 2 ) 


with 


*? = *? = o, 


c _ 1 + a A 
1 1-0Z-V 


*c _ L 1 + a r ^ 4 ‘o -4 « 1 


Zo - zo 


*? = 


^ - Zo 


*2 — *T j 

(z - zo) 3 Z - Z# 


= a ~P \ A . 2 ( z ~ z °) ) 

1 + {} Z — zo (z — zo) 2 ’ 
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(3) 


*? = 


- z&i' + 


*±l A 

1 — j3 z — zo 


A . Azo-zo) 
z-z 0 ^ A ( S ~z 9 r 


In (3) the Dundnrs constants a and j3 are defined by 


+ 1) - + 1) 
” Mi(*a + 1) + Ma(*i + 1) ’ 


g „ p \(*2 - 1 ) - - 1 ) 

1 mi(*» + 1) + + 1) ’ 


(4) 


and 


A = 


_C (1^4 J )j 

27r* (1 — a) ’ 


2/ii(l-ar) 
-(*1 + 1)(1 -/*»)’ 


(5) 


where 5 = 6, + >'6 V with 5, and b v being the x and y components of the Burgers vector. 
We next compute the stresses doe to a dislocation 6i = &i» + ib i„ at the interface (zq — ► 0) 
and doe to a dislocation b? = b-zs located at * = 0, y = ( < 0. The first dislocation 
corresponds to the Green’s function for the interface crack, while the second corresponds 
to the fundamental solntion for the vertical crack. For zq — » 0, the stresses become 


»<>,) + i»m = [|i - ,^cs(*)]4, 


for y = 0, 


( 8 ) 


- £< i + 




( 7 ) 
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for y > 0, and 


•» + *ii> = 7 (p + - f (1 - 0 (^ 






( 8 ) 


with 


1 2y ^ «, a\ V ;tfi (i/j\ * fi o\ 2®y^ i 


P(y,x) = (1- p)?L - (3 - 2/J)i - <1(1 - 2/J)-j - (1 - (>)-£-] 


(9) 


for y < 0, where S(x) is the Dirac delta function. Here, some care must be taken daring 
performing the limit [8]. Meanwhile, for dislocation 63, it can be shown that 

„(!) + *(D = c Mv-Q 2 _ £ + 

a *V T - l r 4 r 2 * ^ -4 

t rj r x 




W - g - r - + 

*1 r l r l 


( 10 ) 


for y > 0, 


<,( 1 ) + ^( 1 ) « 2 {(1 _ fl3)( 2 k~ A 2 * _ £) + (a + 8 2 )( 2 ^ + /fe 

W *(l + a) U PA rf rf ;tl TP A rj 


-4) - (1 + a)/? 2(C y )3g + (a - /?)(! - 
r 2 r 2 ~2 
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for y < 0 , where 


. UKC+ri** ) + 2 (a- 0(1 - - ^))»». 

r 2 r 2 “j 


+* 


. C 


*(1 + o) 


id - - ^ 7 ^) -(«+/»*) 

r f r i 


x( ^.c^ )+(1+a) ,M 


«(C + ») J 


- 2 (« - 0(1 - 0 ci ) - (« - «(1 - 

'2 ~2 


. Mft+y}* ) + J(a _ w _ 0c>(«<L+«!> _ 

r 2 r 2 r 2 


( 11 ) 


r? = x 2 + (y - C) S , r| = x 2 + (y + C) 3 - 


( 12 ) 


In particular, when y = 0 


<r^(x,0) +wr^>(x,0) = ^#(x,C)*j 


( 13 ) 


with 


£(x ,0 = 


*+»C 

*1 


-2(l-/?)x< 


C + « 


(14) 


where r 2 = x 2 + C 3 ? and when y < 0 


( 15 ) 
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where 


F{y,0 = 


(a+P) 1 2(a — ft) C 

(l-/P)y + C 1+ft (y + CP 


4(a -ft) C* 

* + £ (y + C) 3 * 

(16) 

The solution to problem (2) is obtained by using Fourier transformation techniques. 
Without loss of generality, consider a bimaterial strip loaded with surface tractions which 
are symmetric about the y-axis. The Fourier transform of displacement is as follows: 


Ui(£,y) - Ui(x,y)smZxdx 
Vi(£, y) = jf° Vi(x, y)ct*Zxdx 


(17) 


The general solutions for Ui and Vi are [9]: 


Ui({, y) = (A*h 2 + Aay)e~* r + (A*h 2 + A+y)^’ 


V*(t,y) = + (m/Z + yJ-AaJe - ** + [-i4aAa 


+(*»/£ - v)Au}& 


(18) 


with the unknown coefficients determined by continuity and boundary conditions along the 
interface: 
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An + Au — An - A-n = 0 


A u + 


-rr-Au — -Au + 

C »3 


Kl Ax 4 — An — 77 — A 22 


^3 


to 


+A U - 


*3 

toi 


a m = o 


+ 2(1 — pi)Ai 3 + tkjAxa — 2(1 — vi)Ai4 — — 

Mi 


x[f hiAn + 2(1 - J'ajAaa + to -^23 — 2(1 - p^Aat] = 0 


^^3 -All + (1 — 2v\)A\a — ^ 2 -Ai 3 + (1 — 2j*l).Ai4 — — 


x[f&2-An + (1 — 2 wj)A 22 — ZhiAn + (1 — 2v2)An] = 0 


— [f(Aito + A 12 & 1 ) + 2(1 — t'ljA.uJe * hl + [-((A 13 A 1 


+Ai4&x) + 2(1 - vijAuJe** 1 = A(0 


— [f(AnAa + A12&1) + (1 — 2Kx)Ai3]c~^ fcl + [{(Aito + Au^i) 


-(1 - 2^)A 14 lc fA ‘ = MO 


— [f^a(-Aai - Aaa) + 2(1 — Vi)An\^ h ' 1 + [— to(-^33 — -A 24 ) 


+2(1 - i^)Aa 4 ]e-« A » = MO 
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-lto(4u - 4 b) + (1 - 2*a)A B }e**« + [fh 3 (An - A 34 ) 


-(1 - 2v 2 )A 24 ]e-* h ' = A(0 


( 19 ) 


where the first four equations represent the continuity of u, v, tr n and along the 
interface, respectively, the last four equations specify the boundary conditions at y — hi 
and y = —h 3 . The functions fi(0 and MO axe related to the Fourier transforms of the 
loadings on the top surface of the upper strip, while MO and MO axe related to the 
Fourier transforms of the loadings on the bottom surface of the lower strip. To eliminate 
the unwanted surface traction along the boundaries due to each dislocation, for example b z , 
in problem (1), functions MO, MO, MO “<1 MO axe set to be 


MO — f Q *$?(*»** ) co *Z xdx 

MO = ~ £° *$j?(*»*t)«»£*** 

MO = “ 2^7 J *$(*, -h 3 )co*{xdx 


MO = ~ 2 ““ f o — ha) sin dr. 


( 20 ) 


It can be shown that equations (19) are also applicable to the anti-symmetric problem. The 
MO's have been evaluated for each case and axe given in the Appendix. 

After the A,y’s are solved, the stress components due to problem (2), which are denoted 
by superscripts “(2)”, are readily obtained by inverse Fourier transformation. For the 
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symmetric problem, e. g. (or dislocation b a . 


9 ay( x -i 0) = j£°° K^aC^is - ^n) - (1 - 2 j/i)(Al 2 -f Ah)] 


x sm£xd£ 

(*, 0) s* ~~ jf°° [~C&a(An + Au) •+• 2(1 - »*i)(Ai4 - Au)] 


x cos £z d£ 


( 21 ) 


and 

^am (*»y) ® ' jT* {[£(-^21^2 + ^22V) - 2V2A3a]e~^ y 

+[£(•^23^3 + Any) + ZpjAmM*} 008 £ x <^£ 

(22) 

where y < 0. For the anti-symmetric problem, e. g. for dislocation 6„ the above relationship 
is still valid if sin£x is replaced by cos£z and cos(x is replaced by — sin£x. 

For convenience, in subsequent discussions introduce the following definitions for the 
s tr ess e s produced by problem (2): 


*9(*» 0) + *<r$(x, 0) 9 Gi(x)bu + Gi(x)bi v + Gs(x; C)*3» 

•£?(*> 9 ) a H\(v; *)&i* + H 2 (y, x)b lv + H z (y, C, x)bu . 

( 23 ) 
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We wish to point out that the fundamental solutions derived above con be used to solve 
a large class of problems. The purpose of this paper is to present the methodology. Thus 
we will consider only three types of loadings. These include three and four point bending 
(Fig.la), constant pressure along the crack surfaces, and a temperature change of the com- 
posite strip (Fig.ld). Consider first a perfectly bonded bimaterial beam subjected to four- 
point bending. The stresses in the strip due to this loading condition can be calculated using 
composite beam theory if the strip is long enough. However, since we may be interested, in 
the future, in analyzing relatively short beams, an elasticity solution is obtained in the same 
manner as for problem (2) of the dislocation solution. The applied loads are represented 
by Dirac delta functions at the load points. That is <r„(x, hi) = -p[6(x — d) + S(x + d)]/2 
and <r w (x, -hi) = — p[S(x - e) + 8(x + e)]/2. The functions j\, fa, fz and f\ for this 
case are listed in the Appendix. To evaluate the stresses which arise from this loading, it 
may seem natural to proceed in the some exact manner as for the dislocation solutions. 
However, tremendous care must be taken when evaluating the stresses along the interface 
and along the line * = 0 due to this loading condition, because the integrals in the inverse 
transformations converge very slowly. To improve the convergence the integrands are first 
evaluated at £ = oo. The dominant portion of these limits, which correspond to half-space 
solutions, are subtracted from and added to the integrands, and the added integrals axe 
evaluated in closed form. The same procedure is also applied to the dislocation bi when it 
is close to the bottom surface. For brevity the details of the procedure are not given here, 
but can be recovered in similar analyses presented in [10,11]. 

Denote the stresses produced by the four point bending as 


0) + *«**(*, 0) = pQ(x)/hi, 

*••(0, y) s pR(y)/ hi (24) 

respectively, where p is the force per unit thickness. These will be used in section 3 to set 
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up the integral equations. 

The next loading condition consists of a uniform change of temperature AT of the 
composite strip. The thermal stresses in the strip are calculated using beam theory [12]. 
In particular, the stress in layer 2 along the line x — 0 is given in terms of the thermal 
expansion coefficients ati and a-i by 


<r„ = -( 0l -a 1 )Ar£3(l + |X)/( 1 + ^) (25) 

with 


E* = 


E 

1 _„ 2 * 


3 Integral Equations 

Several problem configurations are depicted in Figure 1. The coordinate system is chosen 
such that the x axis lies along the interface of the bimaterial strip, and the y axis along the 
vertical crack so that layer 1 is at y > 0 and layer 2 is at y < 0. The integral equations for 
the four point bending problem shown in Fig.la are set up first. The stresses along the lines 
y s= 0 and * = 0 (y £ 0) axe given lay the summation of contributions from problem (1), 
problem (2) and the four point bend load (24). Replacing the dislocations by a distribution 
of dislocations enables ns to satisfy the traction boundary conditions along the crack. This 
procedure leads to the following set of coupled singular integral equations: 

A f~t^ dt ° ~ ii3 * Bi ® + Ax * i( *’ *° 

+ J ^Ki(t,to)S\(to)dto + j ^ Ki(t, *o)-B2(<o) d£o = pi(t) 

A *0 + r K 4 (t, to)B 2 (te) dto + f Ks(t, *,)*, (* 0 ) d*o 
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( 26 ) 


+ £ J T.(i, to)Bi(to) = M (t) 


together with the conjugate of the first equation in (26) where the dislocation densities are 
given by Bi(t) = Bi K (at) + iB\ v {at) and B^{t) — Bufaii — l)/2) with 


Bl.(x) = J;K*,0+) - .(*,0-)] 

Bi»(x) = -^(.(*,0+) - r(x,0-)) 

B>.(*) = l-WO+.r) - «(0-,j)l (27) 

p °v 

*i(t,to) = g[Gi(o(< - to)) - <<%(«(* - *o))l 
& 2 (Mo) = ^\GMt - to)) + iC? 2 (a(t - 1 0 ))] 

*s(t,to) = ^[£(«t, + £<? 3 (at; ^S&zH)] 

»• f \ _ ^ 2 rp /^ 2 (* ~ 1 ) ^(*0 l )\ *"(1 + a ) 

fir 4 (t,t 0 ) - y[F( - , 2 ) c{i-p) 

»>(«, * 0 ) = -*») + -*o) - »*(»; -*•)]} 

Ka(t,to) = 


14 



pi(t) = -*Q{at) 


M*) = 


1 + or 
l-P 1 


tR( 


h](t ~ 1 ) \ 
2 ' 


(28) 


If the vertical crack is not present, By is set equal zero and only the first equation of (26) 
is enforced. If the interface crack is not present, on the other hand, the first of equations (26) 
is not enforced and B\ is set equal to zero in the second equation. For the edge crack whose 
length c < ky, the dislocation density is instead represented as By(t) = B<u(c{t — l)/2), 
and in the expression of the kernel K\ given by (28) hy needs to be replaced by c. Other 
loading conditions can be treated by modifying functions p\ and py. 

The integral equations can be solved numerically by representing each dislocation density 
in terms of a regular function and a characteristic function with the proper singularities at 
the end paints. Thus, let 


*(*)- 


MO 


(1 - t)7(l + tp-r’ 


By(t) = 

w (1 - t)Vi + 1 


(29) 


where fa (t) and fa(t) are regular continuous functions which are approximated as piecewise 
quadratic [13], and 
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€ 



IzI 

l+0‘ 


(30) 


The exponent A is taken as 0.5 for the T-crack (Fig.la) and for an edge crack whose tip does 
not touch the interface (Fig.lc, Fig.ld). If the tip of the edge crack touches the interface, 
A is determined by the following characteristic equation [14,15]: 
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( 31 ) 






Using the method developed by Miller and Keer [13], integral equations (26) can be reduced 
to a set of algebraic equations. 

The additional conditions to be satisfied are as follows. For the T-crack 



(32) 


fc(±l) = 0 (33) 

j -Bia(t) dt ^ o. (34) 


Equation (32) is the closure condition for the horizontal crack, (33) represents the condition 
that the stress singularities at the tips of the vertical notch are less than square root, and (34) 
implies the bluntness at x = 0. As will be discussed in the next section, the representation 
of the dislocation density for the T-crack is not rigorous. However, the error introduced is 
acceptable for the configurations considered in this paper. 

For a single interface crack 



dt = 0, 


(35) 


and for a single edge crack 


*»(-!)= 0 


(36) 
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J l Butt) it £ 0. 


(37) 


4 Numerical Results and Stress Intensity Factor Analysis 

The stress intensity factor K of the interface crack (Fig.lb) and the T-crack (Fig.la) axe 
defined by [16] 


K = Jim - a)(x - a)- fc [«r w (*, 0) + ur mv (z , 0)]} (38) 


It can be shown that in terms of which is determined by equation (26) with (28) and 
(29), this complex stress intensity factor can be expressed as 


K = gyVl - 


(39) 


for the four point bending load. For an edge crack whose tip does not tonch the interface 


K = , [Jl-xih-i + y - c)a mm (0, y )], 

<e-A,) * 


( 40 ) 


and similarly in terms of 0j(f) 




(41) 


for the fonr point bending. For an edge crack whose tip touches the inter&ce 
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( 42 ) 


K = lim[v / 2Ty' x <r«.(0,y)], 

If— *0 

and 

(43) 

For the other loading cases, the factor p/hj in (39) (41) and (43) is replaced correspondingly. 

The energy release rate of the interface cracks in the considered plane strain problems is 
given by [17] 

<? = [(!- Ki)//n + (1 - n)/(ij]iry/4 co*h 5 (ir<). 

(44) 

In order to check the lengthy algebra and the numerical scheme, two configurations 
which have previously been analyzed are considered first: a pressurized crack between 
two bonded dissimilar layers and a pressurized vertical edge crack in the bottom layer 
of a composite beam. Fig.2 shows the dimensionless stress intensity factors -STyWy/a = 
lim , ,„[(« - x) 1 -'T(a + x) 7 (<r w + i<r av )\J<T oo versus 0.1 < hi /2a < 4 for the interface crack 
with h 2 /hi = 3, v\ = vj = 0.3, = 3 and m/m = 10, respectively. The results are 

plotted in terms of this definition of stress intensity factor for comparison with reference 
[18]. For hi /2a > 0.5 the results agree with those presented in [18] (no results were given 
there for hi /2a < 0.5). As expected, when the crack becomes shorter, the stress intensity 
factor converges to 1 + 2*e which is the solution for an interface crack between two bonded 
half-spaces. 

For the edge crack problem a homogeneous beam and four ceramic composite strips are 
considered. They are: strip 1 Ti/AlaOs with v\ — 0.322, — 0.207 and n 2 / Pi — 4.129, 

strip 2 Ni/MgO with ^ = 0.314, v* = 0.175 and /i 2 //*i = 1.588, strip 3 MgO/Ni and strip 4 
AbOs/Ti [16]. Fig.3 and Fig. 4 show the dimensionless stress intensity factors as functions 
of crack length and relative thickness of the layers for the four point bending case considered 
in [3,4] where d = 5 h 2 , e = 8 jha and I = 3|h 2 . The quantity oy.. T is the axial stress at the 
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lower surface of the composite beam calculated using composite beam theory. It is given by 


with 


_ Mm Hi 

(I i + mli) 


M2(l - 
0l(l - V*)’ 


6 = ki/hi, 


(45) 


H l = k l +hi- H 2 , H 2 = *^ ( ™ + 26 + 6 2 \ 

h = - y) 2 , ^ + W - y) 2 - 

(46) 

where M is the bending moment; for the four point bending shown in Fig.l M — pl/2. It is 
observed that as the crack tip approaches the interface, the stress intensity factor increases 
drastically if the edge crack is in the stiff er layer. In the limit, the stress intensity factor 
approaches infinity for this case, since it can be shown that as the crack tip hits the interface 
the stress singularity is greater than square root. When the crack is in the softer layer, the 
stress intensity factor goes to zero as the tip approaches the interface, since the singularity 
in this case is less than square root. When the crack is very short, namely c/h 2 -* 0, the 
dimensionless stress intensity factor approaches 1.586, which corresponds to the solution 
of an edge crack in a half-space under uniform tension. When the thickness of the upper 
layer becomes very thin, the stress intensity factors for the bimaterial strips approach those 
for a homogeneous beam. It is interesting to note that for these values of mismatch, the 
dimensionless stress intensity factors are insensitive to c/k j for c/h 2 < 0.5 and insensitive to 
k\Jhi for k\/ki > 1. Results were also obtained for three point bending. For this geometry 
the nondimensional results were found to be almost exactly the same as for the four point 
bending. 

Fig.5 shows the stress intensity factors for the edge crack whose tip touches the interface. 
In this case the singularity depends on the elastic mismatch of the two layers and factor A 
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is calculated from (31). For the combination considered here v\ = vi = 0.3, Pi/ Hi = 3 and 
Hi/ Hi = 1/3, and singularities A = 0.6205 and A = 0.4005, respectively. Notice that, since 
the stress singularities are different, it is meaningless to directly compare the values of the 
dimensionless stress intensity factors given in Fig.4 and Fig.5. It should be noted that these 
results do not compare at all with those in [5,6]. 

Snppose that a bimaterial beam is subjected to a concentrated force p along its neutral 
axis as shown in Fig.lf where H\ and Hi are given by (46). The longitudial stress across 
the lower layer is 


<r(2) = 1 

hj( 1 + SS) 


(47) 


where £ = pi(l — vi)I\hi{ 1 — t'l)]. If the beam is subjected to the remote bending, based 
on the composite beam theory, the stress on the lower layer is 


tr t — 0 nux 


Hi — hi — jf 
Si 


(48) 


where <r mu and Hi are given in (45) and (46). Notice that two Dundurs parameters a and 
0 can cover the elastic moduli dependence of a two-dimensional bimaterial system [19]. In 
Thble 1 and Thble 2 the stress intensity factors are given in terms of a and 0 for an edge 
crack with c/hi — 0.5 under the above two loading systems. Table 3 shows the singularity A 
of the edge crack whose tip touches the interface. The corresponding stress intensity factors 
are given for this class of edge cracks. 

Fig.6 and Fig.7 show the stress intensity factors of the edge cracks for the thermal loading 
case. The sign of the stress intensity factor changes when the crack length is approximately 
equal to O.Qhi as a result of bending. If additional loads are superimposed, these results 
suggest that crack growth can either be accelerated when (ai — ai)AT < 0 or slowed down 
when («2 — cti)AT > 0 by the introduction of a temperature change. 
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The last example considers the T-crack, which is referred to as the Santa Barbara speci- 
men, shown in Figure la where i = 5 hj, e = 8 jfcj and l = Z\h^. Results were calculated for 
crack lengths which are shorter than those presented in [3,4]. Figure 8-14 show the dimen- 
sionless stress intensity factors Kk^h^/pl , energy release rates G and the phase angles, 
defined as the argument of Kb.*, for strips with vi =s vi — 0.3 and = 1, Mi /^2 = 2.5, 
= 5, /*i/m 2 = 10, respectively. The results obtained in [3,4] for a/hj > 0.5 axe 
superposed on the figures. The agreement is observed to be quite good. It is observed that 
steady state begins at crack lengths equal to h?. Results for a/hi less than 0.1 were not 
calculated because we believe that accurate solutions cannot be obtained for these cases 
using the present formulation, since the integral equations axe coupled. To obtain accurate 
results for small value of a/k j one would need to formulate the problem in terms of a single 
integral equation whose kernel includes an analytic solution for the interaction of the two 
cracks. 

It should be mentioned that for crack lengths greater than 1.5&2 convergent results were 
not obtained as a result of using a distribution of dislocations density which is continuous 
at x = 0, y = 0. Because the vertical crack intersects the interface crack, the slope of 
the horizontal crack is discontinuous at (0,0). The numerical scheme implicitly assumes 
that the slope there is zero. We believe it is this error which leads to poor results for long 
crack lengths, for which the slope at (0,0) is definitely not zero. This problem is being 
addressed by the authors at the present time by introducing a discontinuous distribution of 
dislocations. 

5 Conclusions 

An analytical formulation has been presented which can be used to solve a class of plane 
clast ostatic problems involving cracked bimaterial beams. Results were calculated only for a 
few geometric configurations and symmetric loadings. However, since the analysis relies on 
fundamental solutions for dislocations in bonded strips, other configurations such as inclined 
or curved cracks and/or non-symmetric loadings can be treated with minor modifications. 
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I For the T-crack, when the interface crack is relatively short, the results obtained using the 

singular integral equation approach compare very well with those obtained in [3,4] using 
the finite element method. However, for such configurations, which involve discontinuous 
dislocation densities at the intersection of the two cracks, poor results were obtained for 
relatively long interface cracks. This problem can be treated by modifying the numerical 
scheme to handle discontinuous dislocation densities. 
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Appendix 


For dislocation b a lying at interface (0,0) 


/i(0=£N 9 + ( 1 + «M]e-‘ ,{ 


MO = -^[i - (i + 


MO = ~IP - (i - WMK* 




For dislocation b v lying at inter&ce (0, 0) 

A« ) = ^[i+ (1 +«MK‘ ,( 
MO = -^\P - (i 
MO = ^ [ i+<i -AM*-** 
MO = — l V+(\ -«ArfIe-«. 

For dislocation lying at (0,() with ( < 0 


(A.1) 


(A.2) 
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MO = - 0( + (h + 

MO * -^[1 - (1 +P)(hi - Ot - 2ptf]c-V"-& 

*0 - ~d?h) [(1 - + <*“*** 

Cb 

^,-; a) W + a) - (a - /J)(l - /J)(A, + (X 
-2(a - «(1 - 

- v^)> ( “ + ' sJ > - - w 1 - «(*• - a 

-2(a - 0)(1 - /S)<ArfV<*»-0* 

(A.3) 

For four-point bending 

/l(0=! '4^ COS ^ A(0 -o 

/*(£) = co»fo A(0 = 0. 

(A.4) 
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Table X Stress intensity factors for edge crack under tension along neutral aads (cfk 2 — 
0.5). 


*/<7< 3 Vire/2 

hx/h? 

P 

a 

-0.8 

-0.6 

-0.4 

-0.2 

ESI 

0.2 

0.4 

0.6 

0.8 


-0.4 

3.891 

3.786 









-0.3 

3.892 

3.788 

3.681 

3.571 







-0.2 

3.893 

3.791 

3.685 

3.576 

3.463 

3.345 





-0.1 

3.894 

3.793 

3.688 

3.581 

3.469 

3.353 

3.231 

3.102 


0.1 

0.0 

3.895 

3.795 

3.692 

3.585 

3.474 


3.240 

3.112 

2.967 


0.1 



3.695 


3.480 

3.366 


3.122 

2.979 


0.2 





3.485 

3.373 


3.131 

2.991 


0.3 







3.262 

3.139 

3.001 


0.4 









3.011 


| Jf/cr< 2 V*c/2 1 



a 



mn 

i|BQ| 

0.0 

0.2 

0.4 

0.6 

0.8 


0.4 

2.782 

2.435 









-0.3 

2.828 

2.481 

2.284 

2.146 







-0.2 

2.866 

2.520 

2.322 

2.181 

2.066 

1.964 





-0.1 

2.900 

2.554 

2.354 

2.211 


1.989 

1.886 

1.773 


1 


2.929 

2.584 

2.383 

2.237 

2.117 



1.786 

1.635 


0.1 



2.408 


2.138 


1.917 

1.794 

1.637 


0.2 





2.155 


1.927 

ESI 



0.3 







1.933 

1.801 

IIS1 


0.4 









1.618 || 


1 KloMyJltcjl || 

ht/h* 

P 

a 

-0.8 

-0.6 

-0.4 

-0.2 

Will 

0.2 

0.4 

BO! 

0.8 


-0.4 

2.117 

1.894 









-0.3 

2.137 

1.910 

1.774 

1.678 







-0.2 

2.149 

1.920 

1.783 

1.685 

1.610 

1.549 





-0.1 

2.157 

1.924 

1.786 

1.687 

1.611 

1.549 

1.498 

1.453 


10 

m 

2.158 

1.924 

1.784 

1.684 


1.545 

1.493 

1.449 

1.410 


0.1 



1.778 

1.677 


1.537 

1.485 

1.440 

1.401 


0.2 





1.588 

1.525 

1.472 

1.427 

1-388 


0.3 







1.454 

1.409 

1.370 


0.4 









1.346 
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Table 2 Stress intensity factors for edge crack under pure bending (c/Ii 2 = 0.5). 
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Table 3 Singularity for edge crack whose tip touches the interlace. 


? 1 

0 

.. 1 

-0.8 

HO 

HD 

n 


0.2 

HO 

HO 

HO 

Bfl 

0.7321 

0.6449 




■ 




wn 

0.7394 

0.6507 

0.5949 

0.5567 


mB 




E * 

0.7438 

EE3 

0.5946 

0.5533 

0.5230 





-0.1 

0.7456 

0.6529 

0.5912 

0.5467 

0.5135 


0.4682 

0.4523 


^ ■ 

0.7450 

0.6495 

0.5843 

0.5364 



0.4496 

0.4318 

0.4173 

nfs 



0.5734 

0.5213 

0.4812 

0.4498 

0.4249 

0.4049 

0.3887 

0.2 





0.4551 

0.4196 

0.3913 

0.3687 

0.3504 

0.3 







0.3435 

0.3173 

0.2963 

0.4 


! 







0.2114 
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Table 4 Stress intensity factors for edge crack whose tip touches the interface (tension 
along neutral axis). 
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Table 5 Stress intensity factors for edge crack whose tip touches the interface (pi 




- 0.6 


12.59 20.29 

10.71 17.87 21.27 22.67 

9.280 16.13 19.47 20.99 21.25 20.56 

-0.1 8.165 14.86 18.20 19.87 20J3 19.85 18.56 16.35 

0.0 7.279 13.94 17.36 19.22 19.90 19.61 18.46 16.33 

16.93 19.03 19.98 19.90 18.87 16.74 

20.73 20.94 20.02 17.81 

22.56 20.14 




hi/fh P I -0.8 I -0.6 I -0.4 I -0.2 


2.478 

2.215 2.536 
2.004 2.324 

2.147 2.260 
1.999 2.130 


0.2 I 0.4 


2.182 

2.077 1.855 


2.130 2.116 1.994 1.790 


1.874 2.024 


KflKSl 


1.745 1.487 1.160 

1.727 1.475 1.148 
1.758 1.505 1.166 


*/<rm«V^(V2) 


hi A* P 


1.734 
1.656 
1.596 
1.554 1.485 

1.531 1.474 

1.489 


0.6 

0.8 

1.274 

1.260 

1.191 

1.263 

1.194 

1.292 

1.220 

1.372 

1.293 


1.497 
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Fig. 8 Stress intensity factors for a T-crack under four point bending {ji\Jp Mi = 1) 
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Fig. 9 Stress intensity factors for & T-crack under four point bending (M 1 /M 2 = 2.5). 
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Fig. 10 Stress intensity factors for a T-crack under four point bending fo/p a = 5 ). 
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fig. 11 Stress intensity factors for a T -crack under four point bending (pi/pi = 10). 
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Fig. 12 Effect of relative thickness of the strips on the stress intensity factors of a j 
T-crack under four point bending. 
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Kg. 13 Energy release rates of a T-crack under four point bending. 
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